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The inner crust of neutron stars consists of nuclei of various shapes immersed in a 
neutron gas and stabilized by the Coulomb interaction in the form of a crystal lattice. 

The scattering of neutrons on nuclear inhomegeneities leads to the quantum correction 
to the total energy of the system. This correction resembles the Casimir energy and 
turns out to have a large influence on the structure of the crust. 

The average density of a typical neutron star is by 3 — 4 times larger than 
the density inside an atomic nucleus. Even though, in the neutron star crust the 
nucleon density is relatively low and do not exceed 0.1 fm“^. The width of this 
layer reaches almost 10% of the star radius and its structure is quite complex. The 
outer part of the crust consists of nuclei forming a lattice immersed in an electron 
gas which become relativistic already at about 10® g cm“®. Deeper in the star, 
in the inner crust, due to high density and pressure, the last occupied neutron 
levels in nuclei are no longer bound, and a neutron gas is formed. The structure of 
this region can thus be viewed as an inhomogeneous neutron matter, where nuclei 
play the role of impurities. The electrons at these densities are ultrarelativistic 
particles and therefore the electron-nucleus correlations, which are responsible for 
the complications of electronic structure in atoms and solids become unimportant. 
Eventually, at the bottom of the crust, the uniform nuclear matter is formed. 

The structure of the outer parts of neutron stars are important for understand¬ 
ing of several observational issues. Namely, one expects that such phenomena as: 
thermal X-ray emission from the stellar surface. X-ray burst sources, or the sudden 
speed-ups in the rotation rate of some neutron stars may be related to its crustal 
properties. It provides thus a motivation for theoretical studies, and indeed the in- 


vestigation of the crv^t structure has been the subject of a considerable theoretical 
effort (see e.g. III’BBu’B and reference therein). In particular, it has been found that 
in the inner crust the interplay between the Coulomb and the surface energies leads 
to the appearance of various nuclear phases, characterized by different shapes. An 
agreement has been reached on the existence of five phases which are formed in the 
region where the nucleon density varies from 0.03 to 0.1 fm“^. Most theoretical 
calculations predict a chain of phase transitions as the density increases: spherical 
nuclei —> rods (“spaghetti” phase) —> slabs (“lasagna” phase) ^ tubes —> bubbles 
—> uniform matter. However in all to date approaches shell effects have been either 
completely neglected, as in the models based on a liquid drop formula or Thomas- 
Fermi models, or have been taken into account only for bound nucleons. Namely, 
the total energy density of the system can be expressed schematically in the form: 


E — Eyoi + E 

surf Eqq^i + Eshell j 


( 1 ) 


where the volume term Ey^ describes bulk properties of the neutron-proton-electron 
(npe) matter, Esurf denotes the surface energy of nuclei, and Ecoui takes into 
account the Coulomb interaction between protons and electrons. The last term 
denotes the shell energy and was usually omitted in the previous approaches. It 
represents the quantum correction to the total energy and has two origins. The 
first one is associated with bound nucleons. Since the nuclei form a lattice the 
single-particle bound states form narrow bands. Their distribution is not a smooth 
function of energy and therefore gives rise to an energy correction which is termed 
shell energy. The change of the total energy can be computed quite accux^ely 
using the shell correction method, once the single-particle spectrum is known B’0. It 
was shown however that this correction play a marginal role and does not influence 
the phase transition pattern 0. The second effect is related to unbound neutrons 
which may scatter on inhomogeneities and form resonant states. The situation 
is somewhat similar tc) tl^ Casimir effect in quantum field theory and condensed 
matter (see e.g. Refs. 00113 and references therein), where the fluctuation induced 
interaction leads to an additional energy correction to the total energy of a system. 
In the same way, the scattering of unbound neutrons leads to an effective interaction 
between nuclei immersed in the neutron environment. This effect can be termed 
lie Casimir effect and has been studied for simple geometries in Refs. 


The shell energy for unbound neutrons can be calculated from the formula: 
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In the above equation gfo(e, 1 ) stands for the Thomas-Fermi or liquid drop density 
of states related to the neutron gas and g{e, 1 ) is the true quantum density of states 
in the presence of inhomogeneities. An ensemble of parameters 1 describes these 
objects and their relative geometrical arrangement. The parameters: /i and go are 



determined from the requirement that the system has a given average density of 
unbound neutrons: 
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Since in infinite matter the presence of various inhomogeneities does not lead to 
the formation of discrete levels, the effects we shall consider here arise from the 
scattering states, which is in complete analogy with the procedure for computing 
the Casimir energy. The only difference comes from the fact that the integration 
in the formula is performed up to the Fermi energy instead of infinity, as in 
the expression for the Casimir energy. This fact is responsible for an oscillatory 
behavior of the shell energy as a function of the distance between nuclei. 

If the impurities were small compared to the neutron Fermi wavelength then the 
generated interaction would originate mainly from the s-wave scattering giving rise 
to the Rudepnan-Kittel correction which is small compared to the Fermi energy of 
neutrons 1130. However the s-wave scattering limit is not valid in this case since 
the neutron Fermi momentum is of the order of Ifm and thus kpR > 1 where 
i? is a size of an inhomogeneity. Consequently the contribution of higher partial 
waves to the scattering process cannot be neglected, giving rise to much larger 
rection, as for large objects more of the incident wave will be reflected 


The determination of ^e shell energy requires the calculation of the scatter¬ 
ing matrix of the system til. It is in principle a difficult task for arbitrary shapes 
and arbitrary mutual geometrical arrangements of inhomogeneities. It is possible 
however to obtain an approximate analytic result under simplifying assumptions 
concerning the nuclear shape. Namely, let us consider here only spherical, cylin¬ 
drical and planar nuclei. Moreover since the amplitude of the wave function in the 
semiclassical limit is proportional to the inverse square root of the local momen¬ 
tum, the single particle wave functions for the unbound states will have a small 
amplitude over the deep well and thus the probability to find an unbound neutron 
inside the nuclei is reduced as compared to the surrounding space? Hence one can 
approximately replace the nuclear potential by an effective repulsive potential of 
roughly the same shape Ii3. Then the shell energy can be easily determined in the 
semiclassical approximation based on the Gutzwiller trace formula. Namely, the 
leading shell energy contribution to the total energy, for two obst^^being either 
spherical, rod-like or slab-like nuclei located at a distance d reads ll3’E2l: 
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where m„ is the neutron mass and i = 0,1,2 for two spherical, cylindrical and 
planar obstacles, respectively (we assume that rods and slabs are parallel to each 

^ There are of course a number of “resonant” delocalized states, whose amplitudes behave in an 
opposite manner. However, the number of such “resonant” states is small and brings only small 
corrections. 







Fig. 1. The total energy density (j^) of the npe matter as a function of the proton quadrupole 
moment Qp = Q 20 ’ calculated in the HF approximation at the constant nuclear density p, and 
proton fraction Z/A, where Z and A are the proton and nucleon numbers per cell, respectively. 
The lattice constant is equal to 26 fm (left subfigure) and 20.8 fm (right subfigure). The stable 
nuclear configurations are indicated by arrows. 


other). In the above equation L defines the length of the obstacle, and R is its 
radius (in the case of a slab it is defined as half of its width). 

Clearly, the interaction induced by shell effects depends on the nuclear shape and 
is a sensitive function of both the neutron density and the geometry of the mutual 
arrangement of nuclei. The interaction between many nuclei of various shapes may 
look, at first glance, quite complicated since three-, four-, and other many-body 
terms will appear as a result of multiple neutron scattering on inhomogeneities. 
One can show however that many-body terms are quto small and give merely a 
small corrections to the dominant pairwise interaction 

The magnitude of the quantum corrections associated with the expression (^) is 
small as compared to the liquid-drop terms in eq. (|^. What is important, however, 
is that the energy differences between various nuclear phases in the crust are also 
very small (of the order of a few keV/fm^), since the liquid drop terms almost 
cancel. One can see it in the figure 1, where the energy density of npe matter has 
been plotted as a function of nuclear deformation. The minima visible in the fig¬ 
ure represent the stable nuclear phases corresponding to different shapes or lattice 
geometries. The results which are presented were obtained through the minimiza- 




tion of the total energy functional (|^ for the npe matter. Namely, we solved the 
Hartree-Fock equations for a fixed density of the npe matter with contributions in 
the energy functional coming from the nuclei, the neutron gas, and the electrons. 
In such a calculation, the liquid drop and the shell energy parts of the energy are 
automatically and self-consistently included. Since we solve the problem in a cubic 
box with three symmetry planes, several lattice geometries can be generated like 
e.g. simple cubi<^a.stal (see), face centered crystal (fee), or body centered crystal 
(bcc) (see Refs. 00 for details). 

Since the liquid drop energy for different phases is almost the same thus the shell 
energy plays a crucial role for determining the most energetically favored structure. 
In particular, it can be shown that the relative energies of different phases fluctuate 
rapidly as a function of the total density and these osc ill ati ons can be attributed to 
the shell effects associated with unbound neutrons Consequently, we may 

conclude that the structure of the crust may be quite complicated since the shell 
effects associated with unbound neutrons may easily reverse the phase transition 
order predicted by the liquid drop based approaches. Moreover the number of phase 
transitions may increase since the same phase may appear for various density ranges. 
It is also likely that the system will favor distorted lattices, or lattices with defects 
which decrease the shell energy. 

Hence our results suggest that the purely quantum effect related to the fluc¬ 
tuation induced interaction between nuclei, and which thus may be termed the 
fermionic Casimir effect, is responsible for determining the structure of the inner 
crust of neutron stars. 
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